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ABSTRACT 



Many chondrules are found to be surrounded by a fine grained rim. Supposedly, these rims were acquired from dust accreted to 
the chondrules in a protoplanetary disk. In numerical simulations we study the heat transfer in illuminated bare and dust mantled 
chondrules. The calculations consider the chondrule size, the dust mantle size, and the thermal conductivities of both components 
as parameters. We calculate the photophoretic force and compare the numerical results to analytical approximations. We give an 
expression to quantify the photophoretic force on a spherical particle in the free molecular regime to better than 2%. We describe 
the influence of a dust mantle on the photophoretic strength by an effective thermal conductivity of the core-mantle particle. The 
effective thermal conductivity significantly depends on the size ratio between mantle and chondrule but not on the absolute sizes. It 
also strongly depends on the thermal conductivity of the mantle with minor influence of the thermal conductivity of the chondrule. 
The size ratio between rim and chondrule in meteorites is found to vary systematically with overall size by other authors. Based on 
this, our calculations show that a photophoretic size sorting can occur for dust mantled chondrules in optically thin disks or at the 
evolving inner edge of the solar nebula. 
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1. Introduction 

The work behind this paper is based on two observations con- 
cerning the mm-size particles, called chondrules, which are a 
major component in primitive meteorites. 

- Chondrules often have specific but different size in different 
meteo rites dKuebler et all Il999b IScott et all 1 1996b iHughesI 
[l978h . 

- Chondrules are oft en surrounded by fine grained rims 
(lMetzleretal.lll99ll) . 

The origin of these facts is still debated. The rims might orig- 
inate in the solar nebula. As bare chondrules move through a 
dusty disk they would collect th e dust. This dust is compressed to 
the observed rims, eventually dBland et al.l 201 1[ Metzler et al.l 
[1991 ICrmeTlt7Dl2008l iMorfUl et al.lll998b ICarballidoll201 ll) . 
From an experimentalists point of view this is very plausible 
as it is almost impossible to prevent that a chondrule in a dust 
rich environment will enshroud itself in a mantle of small grains 
(Bl um & Wurml 120081) . Specific experiments to study accretion 
of dust rims and furth er evolution have recently been carried out 
by iBeitz et all d201 lb . As far as the size sorting is concerned 
also a number of different scenarios are sketched in the litera- 
ture. The local reservoirs of matter in the solar nebula, which 
created chondrules and which fed a parent body, might just 
have produced different size distributions in different locations. 
More often though it is discussed that certain mechanisms might 
have acted on a particle reserv oir, which separated chondrules of 
different size (iLiffm an 2005). Aerodynamic aspects are clearly 
tempting here as pa rticles of different size couple differently to 
the embedding gas dCuzzi et al.lll996l) . All the explanations for 
size sorting and rim formation seem to have their advantages 



and shortcomings. We cannot dismiss one or the other as within 
their set of assumptions all might currently be viable scenarios. 
In this paper we assume that the rims set off as highly porous 
dust mantles on chondrules of nebula origin and add just one 
more assumption, namely that the particles are illuminated by 
the young sun. Adding a light source might then naturally lead 
to a (size) sorting due to photophoresis. 

Photophoresis w as first introduced to th e physics of pro- 
toplanetary disks by iKrauss & Wu rm (2005). First attempts to 
explain size sorting and concentration of chon drules by pho- 
tophoresis were made bv lWurm & Kraussl (|2J306j). They showed 
that in an optically thin solar nebula chondrule aggregates are 
sorted due to the size of the constituents in different locations. 
Key factor here is the thermal conductivity of the aggregates 
which determines the temperature gradient upon illumination. 
First experiments t o measure photopho resis on chondrules have 
been carried out bv lWurm et al.l(l2010l) . The measurements sup- 
port the idea that photophoresis acts sufficiently strong on chon- 
drules to provide sorting. The assumption of a directed radia- 
tion source is vital to the problem and has two very clear sides. 
On one side, without a directed radiation source (visible or in- 
frared) photophoresis does not exist. Therefore, photophoresis 
does not exist in the midst of an optical thick protoplanetary 
disk. However, on the other side, given a light source it is also 
inevitable that photophoresis acts as significant force on the 
particles. In the inner parts of an optical thin protoplanetary 
disks it is orders of magnitudes larger than a stars gravity and 
cannot be neglecte d dMoudens et al] 1201 lb IKrauss et al.l 120071: 
[ Mousis et a0 l2007b IWurm & Haackl 120091: iTakeuchi & Kraussl 
l2008HHerrmann & Krivovll2007l) . The assumption that the solar 
nebula was optically thin at certain times and places is supported 
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by observations of protoplanetary disks. There are a number of 
disks observe d now which have large inner clearings, which are 
optically thin (ICalvet et al.ll2002t ISiciha-Aguilar et al.ll2008l) . In 
some of these disks the star is still accreting matter so that the in- 
ner region still must contain gas. Clearly the conditions for pho- 
tophoresis are given here. Also, disks in general have an inner 
edge. At the inner edge of a disk sorting of chondrules is possi- 
ble as some particles might be pushed into the d ark a bit stronger 
than o ther s. This idea has be en suggested by lHaack & Wurml 
d2007l) and lWurm et alJ d2006l) . To evaluate the potential of size 
sorting and photophoretic transport in further detail in the fu- 
ture we carried out numerical heat transfer calculations on bare 
chondrules and dust mantled chondrules as prerequisite here. 
These are used to quantify the photophoretic strengths on par- 
ticles more accurately. The structure of this paper is as follows. 
Section 2 describes our heat transfer simulations, which result in 
temperature fields at the surface of a (dust mantled) chondrule 
that is illuminated. In section 3 the resulting temperature field at 
the surface of a particle is used to calculate photophoretic forces. 
These forces are compared to values, which are calculated from 
common equations found in the literature. To facilitate analytical 
treatment of photophoretic forces in future work we introduce 
correction factors on the order of 20% for homogeneous par- 
ticles of given size and thermal conductivity. These results are 
then used in section 4 to define average thermal conductivities 
for dust mantled chondrules depending on the chondrule size, 
dust mantle size and the thermal conductivities of both com- 
ponents. In section 5 we sketch how chondrules might be size 
sorted due to photophoresis based on our calculations. 
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2. Heat transfer calculations 

Our basic model setup consists of a spherical particle, which is 
made up of a core (bare chondrule) and a dust mantle (rim) (fig. 

Both components have specific thermal conductivities. The 
core, being compact and consisting of silicates, ha s a thermal 
condu ctivity on the order of k core = lW/(mK) (lOpeil et alJ 
120101) . The dust mantle being highly porous has a very low 
thermal conducti vity. We vary kd ust between 0.01 W/(mK) 
and 0.5 W/(m K) dPreslev & Christensenl[T997allb1: iKrause et all 
l201lh . This core-mantle sphere is illuminated from the top (fig. 
Q]), absorbs radiation and heats up until in equilibrium it emits 
as much thermal radiation as it absorbs visible radiation. Due to 
the directed illumination a temperature gradient over the particle 
surface is established. The problem to solve is the heat transfer 
equation 



V • kVT = Q 



(1) 



where k is the local thermal conductivity, T is the local tem- 
perature and Q is a volumetric heat source. As boundary condi- 
tion we have cooling due to thermal emission and heating due 
to illumination. We consider a complete absorption of incident 
light at the surface (i.e. Q = 0). We further assume an emissivity 
for thermal radiation of 1 . For a plane wave light source of light 
flux / * e/ we then get a von Neumann boundary condition 



Fig. 1. Model setup of a chondrule core (high thermal conduc- 
tivity) and a dust mantle (low thermal conductivity). The light 
source is directed downwards. Due to rotational symmetry only 
the temperature distribution along one line from top (180°) to 
the bottom (0 °) has to be considered. 



to solve differential equations (COMSOL v4.1 from COMSOL 
AB, Sweden). In case of no rotation the temperature distribu- 
tion is symmetric around the axis given by e/. The maximum 
temperature T max is at the top surface element facing the light 
source (at 180°). The minimum temperature T m - m is at the op- 
posite side (bottom or °). The environment is assumed to have 
T gas = 293 K. This is far below the equilibrium temperature for 
an illuminated par ticle but reflects experimental settings, e.g. in 
IWurm et"aT1 (120101) . It will also allow us to specify dependencies 
on temperatures more detailed and not just refer to a single av- 
erage gas temperature, i.e. we will not assume one average tem- 
perature for the gaseous environment and the particle. A sample 
of calculations is shown in Fig.[2]for a homogeneous 1.12 mm 
particle (total diameter). To calculate homogeneous particles we 
use the same thermal conductivity for the dust rim and the chon- 
drule (core). 

All areas of validity and restrictions of the model are dis- 
cussed in Sec. [6] 



ytvrn = /e 7 



n _ ^4 _ T 4 



) 



(2) 



where n is normal to the surface pointing outwards and the 
first term on the right-hand side is only considered for the upper 
illuminated half sphere of the particle. The light flux in all calcu- 
lations is fixed at / = 20kW/m 2 . Eq.[T]and[2]are solved numeri- 
cally using COMSOL, which is based on a finite element method 



3. Photophoretic forces 

3.1. general case at low gas pressure 

The photophoretic force depends on the temperature distribu- 
tion over the surface of the chondrule. We consider the case of 
free molecular flow or low gas pressure here, where the mean 
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Fig. 2. Temperature distributions for homogeneous spheres of 
1.12 mm diameter with varying thermal conductivities. 



free path of the gas molecules is much larger than the parti- 
cle size. In that case the interaction between particle and gas 
can be treated as individual collisions of the gas molecules with 
the surface. A gas molecule impinging the surface of the par- 
ticle is either reflected specular or diffusely. In the first case, 
due to the larger mass of the solid particle the absolute momen- 
tum of the gas molecule does not change and does not depend 
on the temperature of the particle. If this first fraction of par- 
ticles, 1 - a, is the same over the surface of the particle, then 
due to symmetry, specular reflected gas molecules on average 
do not transfer momentum to the particle. The second fraction 
of gas molecules, a, which are reflected diffusely, stick to the 
surface for a short while, accommodate to the surface temper- 
ature and the momentum of the reemitted gas molecules dif- 
fers from the momentum of the incoming gas molecules. On a 
hot surface part the gas molecules leave faster than on colder 
parts. Therefore the gas carries a net momentum after the inter- 
actions. The total force on the particle is the integral over the 
total particle surface of the local pressure induced by the inter- 
action at the local t emperature and in the free molecular flow 
regime is given as (Rohatschek & Zulehnerlfl985l IRohatschekl 
1 1 9951: iHidv & Brockll 1 97oT " 



F = -i 
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( 






§ p 


1 + 




do- 



(3) 



where 



T' = T gas + a[T - Tg as ) 



denotes the statistical mean temperature and a is the thermal 
accommodation coefficient describing the fraction of gas parti- 
cles being adsorbed by the solid particle's surface, T is the local 
temperature of the particle surface, T gas and p are the gas tem- 
perature and pressure far away from the particle, respectively. As 
Tgas in our calculations we use 293 K. 

3.2. rotational-symmetric temperature field 

In case of a spherical rotational-symmetric temperature field T = 
T(ff), e.g. if a sphere is illuminated with an intensity / in direction 
e/, without loss of generality let e/ := -e, (see fig.Q]) in order to 



develop T in Legendre Polynomials P„: 

T(9) = Y j A„P„(cos(6)) 



(4) 



with Aj as the / expansion coefficient. iRohatschekl (1 19951) fur- 
ther linearizes the square root to its first order. Due to the 
Legendre Polynomial's orthogonality relation F (see eq. [3} is 
only a function of A \ ; all other Aj do not show up sin ce the corre- 
spond ing integrals yield zero. Therfore, eq.[3]yields (IRohatschekl 
119951) 



n p 
-—a — 

3 T m 



■r 2 Aie z 



(5) 



where r denotes the sphere's radius. 

Since A\ is usually unknown there are three well known ap- 
proximations for this coefficient. The first one is (indication (;) 
refers to the i' h approximation presented) 



A\ l) ^ -AT , 
1 2 



(6) 



with AT being the maximum temperature difference between 
top and bottom AT = r max - r m j n . Thus eq.[5]yields an approxi- 
mated force F^ for a non-rotating illuminated sphere: 



—a r AT . 

6 Tg as 



(7) 



Our first goal here is to calculate the photophoretic force for 
a homogeneous sphere (core and rim share the same thermal 
conductivity k) based on eq. [3] exactly but describe the results 
in a way similar to the simple analytic expression of eq. [7] In 
the following we therefore determine a correction factor from 
simulations. We define as correction 



F 



(8) 



Varying size r and thermal conductivity k we find that in gen- 
eral, the approximation overestimates the photophoretic force. 
The correction factor can be considered a constant of z <l> = 0.7 
for thermal conductivities larger than about k = 0.2 W/mK in- 
dependent of the particle size (fig. [3J. For smaller thermal con- 
ductivities there is a dependency of on both parameters. As 
the thermal conductivity of dust falls in the low range a non- 
constant correction factor has to be considered. The parameters 
studied range from r = 0.1 ... 11 mm, k = 0.01 . . .8 W/mK, 
and a — 0.1 ... 1, the latter is explicitly part of eq. [5] To give an 
analytic equation we tried different simple expressions. One ex- 
pression that fits the correction factor over the whole range with 
an accuracy being better than 2% is 



0.740 - 0.209e 



-0.00435 f 



+ 0.478e" 104ff 



(9) 



It should be noted, that in eq.l9l fT3lfT8l and[T9lr is in units of mm 
and k in units of W/(m K). 

The next step is expressing the temperature difference in 
terms of an illumin ating flux. For a ho mogeneous sphere of ther- 
mal conductivity k IRohatschekl (1 19951) finds 



1 Ir 
2~k 



(10) 



where I is the radiative flux. As IRohatschekl (1 19951) notes, 
this is an approximation for good conductors. With low ther- 
mal conductivity dust mantles on high conductivity chondrules, 
it is not clear that the approximation of eq. [10] is sufficient to 
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r (mm) 




r t2 > 



= 0.590 - 0.0599e 



-0.0359? 



+ 0.449e 



-0.997a 



(13) 



This correction factor (fig. ID can be rather far from unity es- 
pecially for low thermal conductivity particles. Clearly, particles 
with thermal conductivities much lower than lW/mK can no 
longer be regarded as good conductors in the sense of eq. [10] 
Doing so would usually overestimate the photophoretic force 
significantly. This is due to the fact that the approximation of 
eq. [10] neglects the effect of thermal radiation. This aspect can 
be improv ed by adding a t erm to eq.QJj] which accounts for this 
(e.g. cmp. iBeresnev et al.l (119931) ) 



Fig. 3. Ratio z (1) = for a given particle size and thermal con- 
ductivity k. For thermal conductivities larger than k = 0.2 the 
factor is essentially constant z (1) = 0.7 (a = 1, this plot) (resp. 
z (1) = 0.8 for a = 0.8). 



study the details of photophoretic particle transport and sorting. 
However, eq.[3]and eq.QJjJare useful as they provide a rather sim- 
ple analytical treatment of photophoresis including all important 
quantities, i.e. the gas properties p, T gas , the radiative flux / and 
accommodation coefficient a, and the particle properties r, k. To 
keep this simplicity but to be accurate enough to apply the equa- 
tions to problems of photophoretic transport of chondrules we 
chose again a correction term. We extend our correction factor 
including deviations introduced by the approximation of eq. [10] 
and define 



2 * + 4oT 3 

r ' ^ u 1 par 



(14) 



Usually, particles are considered in thermal equilibrium 
with their surroundings, photophoretic temperature gradients are 
small and illumination is not significantly changing the average 
temperature. In that situation the particle temperature T par is ap- 
proximated by T gas . For our calculations and often for experi- 
mental settings this is not the case, i.e. the particle temperature 
can easily be a factor 2 higher than the surrounding gas tem- 
perature. As the radiative cooling is depending on the particle 
temperature it is the particle temperature and not the gas tem- 
perature, which enters in eq. [14] As seen in fig. [2] the particle 
temperature varies strongly over its surface and is still not well 
defined. We use the average black body temperature here and 
define 



with 



7 (2) = J_ 
p(2) 



n p 3 J 
—or— — r - 

O I ens K 



(ID 



T = 

1 nar — 



— + T 4 

4a gas 



(15) 



Based on eq.[T4land[T5lwe finally consider as analytical ap- 
proximation that captures all physical parameters of the system 
but is accurate enough to treat dust mantled chondrules 



(16) 




0.85 



Jlh 



Fig. 4. Ratio z (Z) = j^rj f° r a given particle size and thermal con 
ductivity k. For thermal conductivities larger than 2 the factor is 
essentially constant. The accommodation coefficient a is 1 . 



As before, the approximation (eq.[T2l overestimates the pho- 
tophoretic force following the same pattern as before. In the 
same parameter set (r = 0.1 ... 1 1 mm, k = 0.01 ... 8 W/m K, 
a = 0.1 ... 1) z (2) can be well fitted again to better than 2% by 
the following expression 



Fig. 5. Ratio z ( ' = j^j for a given particle radius and thermal 
conductivity k. For thermal conductivities larger than 0.2 the fac- 
tor can be regarded as constant. The accommodation coefficient 
a is 1. 



and quantify the photophoretic force adding a correction fac- 



tor Z 



(3) 
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F = z P) F (3) . (17) 

Similar in structure, z (3) , can well be approximated in the 
whole range of parameters to better than about 2% by the fol- 
lowing terms 

z (3) = 0.725 - 0.177e-° 00205 s + 0.434 e - 114ff . (18) 

As summary, if the temperature gradient is known the first 
approximation (eq. [7j and [3 might be used. More general and 
used further in this paper is the photophoretic force on a homo- 
geneous, perfectly absorbing spherical particle as given by eq. 
[l6]to[T8]or combined as 

F = (0.380 - 0.0927e-° 00205 J + 0.227 e - U4ff ) * 

* ff JL r 2 / _ (19) 

Tgas + T UY 




Fig. 6. k over ddust and kdust for a 1 mm sized core with a ther- 
mal conductivity of 4 W/mK . The accommodation coefficient 
a is 1; white lines indicate constant values and show multiple 
configurations can lead to the same k. 



4. Dust mantled chondrules 

In view of equation [T9j thermal conductivity of homogeneous 
spheres is a major particle property that determines photophore- 
sis. For dust mantled chondrules there are two distinct differ- 
ent thermal conductivities for the core and the mantle and it is 
not clear, a priori, what determines the photophoretic strength of 
such particles. A core mantle particle can be described by five 
parameters, which we vary within the following range: 

- the core radius, 
r core : 0.2 . . . 1.0 mm, 

- the mantle thickness, 

ddust- 0.1 ... 1 in fractions of r core 

(the total radius then is p := r core (1 + ddust)), 

- the thermal conductivity of the core, 
k core : 1 ...4 W/mK, 

- the thermal conductivity of the mantle, 
W 0.01 ...0.5 W/mK, 

- and the accommodation coefficient, 
a: 0.1... 1. 

For each combination of these parameters a photophoretic 
strength can be calculated. In principle we could argue on 
the photophoretic strength depending on these parameters 
at this point. We note that different configurations y = 
{r core , ddust, k C ore, k-dust) of dust coated chondrules might result 
into the same photophoretic strength, even if the temperature 
field across the mantle is absolutely different. So two configu- 
rations y\ , 72 are related to each other if 

Ti ~ 72 :o F(y{) = F(y 2 ) . (20) 

Nevertheless, it is possible to choose a representative y (p, k) of 
all those configurations obeying eq. [20] with a (total) radius p 
and constant thermal conductivity a:. We call k effective thermal 
conductivity. 

In view of quantifying thermal properties of dust man- 
tled spherical chondrules with a total radius p(r core , ddust) = 
feme (1 + ddust) we use the results of the photophoretic 
strength calculations and eq. [16] and [18] to calculate 
k (r core , ddust, kcore, kdust)- It is of course not injective, so multiple 
configurations can lead to a same k (see fig.©. In the following 
the dependency of k on the different parameters is considered. 
We therefore numerically solved 



F (T, a) = z (3) (p, k, a) F <3) (p, k, a) (21) 

for k, where T — T (r cor e, ddust, k CO re, kdust) is the sphere's sur- 
face temperature (depending on all four parameters for the heat 
equation system). 

Accommodation coefficient: The accommodation coeffi- 
cient is not connected to k (by construction). Numerically, vari- 
ations only influence the result on the 1 % level. This is in agree- 
ment with the idea that thermal conductivity is independent of 
the photophoretic property. 

Core and mantle size: In view of estimating the potential for 
size sorting the most important dependencies to explore are the 
influence of the total size of the chondrule and the influence on 
the ratio between mantle thickness and core radius. We find that 
only for low thermal conductivities of the dust mantle does the 
effective thermal conductivity slightly depend on the total size. 
The extreme for a dust mantle of the lowest thermal conductiv- 
ity studied of kdust = 0.01 W/mK is seen in fig. [7] The thermal 
conductivity of the core is 1 W/mK. All three lines vary very 
little. This shows that the dependency on total size for a given 
mantle to core size ratio is very weak. For given thermal con- 
ductivities of core and mantle only the ratio between the sizes 
of both components determines the effective thermal conductiv- 
ity. This trend can be understood as the smaller the mantle the 
more effective is the high thermal conductivity of the core. In the 
extreme of a vanishing mantle the thermal conductivity has to 
reach the core value, which is consistently confirmed by the cal- 
culations. For thick mantles the low thermal conductivity of the 
mantle should dominate, which is observed in fig. [7] However, 
it should be noted that even for equal size components the effec- 
tive thermal conductivity is still a factor 2 larger than the thermal 
conductivity of the pure dust mantle. 

Thermal conductivity of core and mantle: From fig. [7] it 
is obvious that the thermal conductivity of the mantle is a dom- 
inating factor for the effective thermal conductivity. In fig. [8] the 
dependency on the size ratio between rim and core and the ther- 
mal conductivity of the mantle is shown for two core thermal 
conductivities. The two different plots are based on the extreme 
values used for the thermal conductivity of the core in our calcu- 
lations. The comparison illustrates that the core has a significant 
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Fig. 7. Effective thermal conductivity k over dust mantle to core 
size ratio. The 3 different lines represent 3 different core sizes: 
from bottom to top r core = 0.2,0.5, 1.0 mm The thermal con- 
ductivity of the dust mantle is 0.01 W/mK, the core's thermal 
conductivity is 1 W/m K and the accommodation coefficient a 
is 1. 



influence only for thin mantles and then only if these mantles 
have high thermal conductivities. 



k dust (W/(m K)) 




K(W/(m K)) 



Fig. 8. k over dd us t and kdust for a 1 mm sized core with a ther- 
mal conductivity of 1 W/mK (bottom) and 4 W/mK (top). The 
accommodation coefficient a is 1 . 



Fig. [9] shows the dependency in higher resolution again by 
plotting k over ddust but for two fixed thermal conductivities of 
both the mantle (0.01 W/mK, 0.1 W/mK) and core (1 W/mK, 
4 W/m K). Even for thin dust mantles, if the thermal conductiv- 
ity has a reasonable value attributed to dust, then the effective 
thermal conductivity does not significantly depend on the ther- 
mal conductivity of the core. 

Conclusion section 4: In total we find that the effective ther- 
mal conductivity significantly depends only on three of the four 
parameters, the relative thickness of the dust mantle and the ther- 
mal conductivities of core and mantle but the total size plays 
no major role for thermal conductivity. The thermal conductiv- 
ity of the core is visible in the values of the effective thermal 
conductivity, i.e. by increasing the value of the mantle by a fac- 



1 W/(m K) 
1 W/(mK> 
01 W/(m K) 
01 W/(m K) 




Fig. 9. k for different pairs (k, 
0.01 W/m K, two top lines: kdust 



^, kdust)- Two bottom lines: k dus t 
0.1 W/mK 



tor of two or more. However, the effective thermal conductivity 
is insensitive to variations of the core thermal conductivity in 
the range of reasonable values of kdust- Assuming that the solid 
chondrules (core) have a fixed thermal conductivity on the or- 
der of 1 W/m K, this leaves two variables in the context of pho- 
tophoretic forces on chondrules, the thermal conductivity of the 
mantle and the dust mantle to core size ratio. 

5. Size sorting of chondrules 

The dependency of the effective thermal conductivity on the dif- 
ferent parameters provides various options by which dust man- 
tled chondrules can be sorted due to photophoresis. We only 
give a qualitative account here. Detailed treatments would re- 
quire specific disk models. This would be a problem based on 
the results of this paper but o therwise being beyond i ts scope. 

The sorting mechanism in 1 Wurm & Krau ss (2006) was based 
on a disk density profile which decreases with distance to the 
star. This way, the gas is supported against the stellar gravity 
and rotates slower than with Keplerian speed in order to ro- 
tate on a stable orbit. Small solid particles like chondrules em- 
bedded in the disk would couple to the gas rapidly, would ro- 
tate with the same slower orbital frequency but would not feel 
the s upport by the pressure gradient. Such particles drift in- 
ward dWeidenschiirind["l977l) . The particles are subject to a ra- 
dial force 



F D (a) = 



nRr, T, 



G -*- gas ™par 

fia 



(22) 



then which can be considered as residual gravity. Here a is the 
radial distance to the star, ji is the molar mass of the gas, n 
is the power of the radial distance dependency of the gas den- 
sity (» = 11/4 in th e minimum mass solar nebula, see eq. 1271 
dHavashi et al.ll 19851) 1 Rc is the gas constant and m par is the par- 
ticle mass. This includes gas drag in the way that the gas grain 
coupling time is much shorter than the orbital timescale and the 
pro blem reduces to a one di mensional radial problem. 

IWurm & Kraussl (120061) consider the simple photophoretic 
force estimate given by eq. [12] and calculate the ratio 



(23) 



If this ratio is 1 inward and outward drift balance and chon- 
drules are concentrated. This concentration distance depends on 
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Fig. 10. Principles of sorting dust mantled chondrules by pho- 
tophoresis 



the strength of the photophoretic force. Therefore, chondrules 
with different photophoretic properties - as outlined in this pa- 
per - would concentrate at different locations. 

Particle sorting due to photophoresis might also take place at 
the i nner edge of the disk, which might move over the asteroid 
belt (lHaack & Wurmll2007l) . The idea is sketched in fig. [TO] At 
the inner edge of a disk which is optically thick the sorting might 
work differently than described above. Per definition the optical 
depth decreases further away from the edge. Particles are there- 
fore pushed outward as long as the decreasing photophoretic out- 
ward motion is more effective than any other inward directed ef- 
fect like turbulent diffusion or any other remaining radial inward 
drift and the particles pushed more vigorously by photophoresis 
will move outward into the dark the furthest. Here, the sorting 
distances are strongly disk dependent but the principle is inde- 
pe ndent of the detail s. 

IWurm & Kraussl (120061) considered that particles with ther- 
mal conductivities varying by a factor 2 could easily be concen- 
trated at different positions in a solar nebula almost 1 AU apart 
of each other. A rough linear scaling would imply that a differ- 
ence in effective thermal conductivity on the order of 10% would 
separate particles on the order of 0.1 AU or 15 x 10 6 km. This is 
certainly sufficient to separate particles far enough to be incor- 
porated into different parent bodies, eventually. 

Direct size sorting: As calculated above the effective ther- 
mal conductivity does not change with total size. The pho- 
tophoretic force has an intrinsic dependence on size though as r 3 
(eq. |T9l > with the correction factor z (3) included. For a homoge- 
neous particle the gravitational pull by a star also depends on r 3 
and the ratio of both forces is essentially independent of the to- 
tal size. Without correction factor photophoresis would not sort 
by size directly. For the lowest thermal conductivities there is a 
size dependence in the correction factor which varies the pho- 
tophoretic force by 10% for the smallest and largest chondrule 
(see fig- [5j. In that case direct size sorting of dust mantled chon- 
drules with plausible mantle thickness by photophoresis is pos- 
sible. It might also be considered here that the dust mantle has a 
different porosity as the core, which might vary by about a factor 
of 3. Nevertheless, for a constant size ratio between mantle and 
core a change in total size would not change the average density 
of the particle and this is not important. 

Mantle to core size ratio sorting: It is interesting to note 
that for a constant dust mantle to core size ratio no sort- 



ing due to the effective thermal conductivity would occur. 
How ever, considering the growth of dust mantles by accre- 
tion, ICarballidol d201 ll) shows that the size ratio between man- 
tle and core changes with overall size. This introduces an abso- 
lute size dependence of the effective thermal c onductivities then : 
K(ddust(p))- Values estimated from fig. l9l in ICarballidol d201 ll) 
are e.g. a 0.3 mm core with a 0.2 mm mantle {ddusi = 0.6) but a 
1 mm core only has 0.4 mm mantle (dd ust = 0.4). Corresponding 
effective thermal conductivities from our simulations assum- 
ing the core to have, k core = lW/mK and the rim to have 
k dust = 0.01 W/mK are k = 0.02 W/mK and k = 0.03 W/mK, 
respectively. This is a difference of a factor 1 .5 which is certainly 
sufficient to separate these particles by photophoresis. A thinner 
mantle also corresponds to a slightly higher average mass den- 
sity and gravitational attraction. This works in the same direction 
as decreasing the strength of the photophoretic force. 

Mantle thermal conductivity: The exact thermal conduc- 
tivity of the mantle is unknown as it depends on many param- 
eters. Assuming that the material properties of the dust are the 
same on averag e, porosity of the du st rim is still a free parame- 
ter. The work bv lBland etaT](l201lh suggests that mantle porosi- 
ties are on the order of 0.7 or more. This is consistently on the 
same order of poro sities to be expected from dust experiments 
( Teiser et al. 1201 ll). For this porosity and micrometer size dust 
iRrause et al 1 1201 ll) find thermal conductivities on the order of 
0.01 W/mK and below. We assumed this value for dust mantles 
above. Typically smaller particles collide wi th dust at lower ve- 
loci ties than larger ones . The experiments bv lTeiser et alj d201 ll) 
and iKothe et al.l (l20Tlh show that such differences in collision 
velocities produce a difference in porosity. In that case the dust 
mantle porosity would systematically change (decrease) with 
total size. Typical differences are several percent. Assuming a 
porosity variation e.g. of 0.05 would lead to a difference in 
thermal conductivit y on the order of 1.5 for the dust mantle 
dKrause et alj|201 lb . Varying the dust mantle thermal conduc- 
tivity by this factor, e.g. between 0.01 W/mK and 0.015 W/m K 
for a core with 1 W/mK changes the effective thermal conduc- 
tivity by a similar factor of 1.54 (r core = 1 mm, ddust = 1) or 1.46 
(fcore = 1 mm, ddust = 0.5). This is on the order of the changes 
induced by the different mantle to core size ratios discussed in 
the section before. This provides the basis for a photophoretic 
size sorting. The trend would go in the same direction as the size 
ratio variations, i.e. larger particles have a lower effective ther- 
mal conductivity and would be found closer to the sun. 

Core thermal conductivity: The third parameter influenc- 
ing the effective thermal conductivity is the thermal conductiv- 
ity of the core. As it is plausible that dust mantles have rather 
low thermal conductivity, the dependency on this parameter is 
not strong. As found in the simulations effective thermal con- 
ductivities for thin (20% of core radius) dust mantles vary by 
10%. For more conductive rims of kdust =0.1 W/mK this might 
be somewhat larger but compared to the dependency on the size 
ratio between mantle and core the effect is small. We conclude 
that for core-mantle particles the thermal conductivity of the core 
wil l not be responsible for size sorting. 

IWurm & Kraussl d2006l) estimate the drift velocity of chon- 
drules due to photophoresis at 1 AU which would be on the or- 
der of 0.1 m/s. Separating the particles by 10 x 10 6 km would 
require about 10 7 s or 1 year or 1 orbit then which is fast com- 
pared to evolutionary time scales. If all mass is in chondrules 
and if chondrules do not evolve e.g. by growing to larger bod- 
ies then a maximum mass of separated chondrules can be esti- 
mated by the number of chondrules which would make the disk 
optically thick themselves. Assuming a disk height of 0.1 AU 
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at 1 AU distance to the star the area of the inner edge would 
be 2n x 1 AU x 0. 1 AU. Completely covered with chondrules 
of 1 mm diameter, the number of chondrules is 6.35 x 10 27 and 
assuming these chondrules have a mass density of 3 g/cm 3 one 
chondrule has a mass of 1 .6 mg. The total number then is 10 22 kg 
which is about 10 times the mass of Ceres. Being that large the 
estimate shows that the mechanism provides more than sufficient 
separation rates even if conditions are less perfect than assumed 
for this estimate. 

6. Caveats 

Radiation pressure: Photophoresis is driven by radiation. Other 
movements of particles induced by radiation and more com- 
monly considere d are radiation pressu re and the Poynting- 
Robertson effect dKlacka & Kocifaill2002l) . The difference is that 
radiation pressure is based on the direct momentum transfer of 
photons while photophoresis is based on momentum transfer via 
gas molecules. For chondrule size particles and plausible gas 
densities photophoresis is much larger as seen in the following 
estimate. Radiation pressure can be approximated by 



F, 



■ad 



(24) 



Compared to the classical photophoretic force estimate of eq.[T2l 
which is sufficient for this estimate as photophoretic force it is 



F w IF rad = 



cr p 
6 k T„, 



(25) 



Taking a typical value of k = O.lW/mK, a worst case of a 
small chondrule of r — 0.2 mm and a worst case temperature 
of T gas = 280 K we get F„h nt /F rn j ^ 3 57 p. The minimum mass 
solar nebula by lHavashi et alj d 19851) has an ambient pressure 
of about 1 .4 Pa at 1 AU. Here, photophoresis is a factor of 500 
stronger than radiation pressure. The pressure is 0.04 Pa at 3 AU 
or distances related to the outer asteroid belt. Here, photophore- 
sis is 14 times stronger than radiation pressure using the above 
conservative estimates. Therefore, photophoresis dominates as 
radiative force. 

Rotation: Particle rotation was not considered here though 
it is important in the context of photophoresis. If particles would 
rotate more rapidly than the typical heat transfer time then the 
temperature gradient would not reach equilibrium and be re- 
duced which leads to a decrease in the photophoretic force. We 
did not consider rotations for the following reasons. Any ran- 
dom rotation is damped on a time scale of the gas-grain coupling 
time. In the free molecular flow regime which is relevant for 
m m-particles in prot oplanetary disks the coupling time is given 
as dBlum et al 

m par 1 

t c = e , (26) 

v par Pgas^gas,rh 

with p gas being the gas density, <T par the particle's geometrical 

cross section, and the mean thermal velocity v gas ,t h - -^ 8 ^'°' - 
At a = 1 AU in a minimum mass nebula it is dHavashi et al] 
119851) 

.-11/4 



p gas (a) = 1 -4 ■ 10- 6 ( j^j) kg/m 3 , (27) 
««) = 280 — K, (28) 



p(a) = 3551 p gas (a) T gas (a) — ^ 



(29) 



Then, r c is on the order of 10 3 s for a 1 mm particle. Random 
rotation is therefore decaying rapidly. Random rotations would 
only be important in decreasing the photophoretic strength ef- 
fectively if they would continuously be excited on similar time 
scales . A way to generate particle rotation is inter-particle colli- 
sions. iKrauss et al.l (120071) estimated that the time between col- 
lisions of mm-size particles of mass-density p par in a 'typical' 
disk is given by 

Tcofl = 3 3^— , (30) 

Pgas *par,rel 

which then is about a factor of 10 4 longer than the gas grain cou- 
pling time (Vpar.rei = 1 m/s). Therefore, even assuming strongly 
varying disk conditions, rotation is damped much more rapidly 
and the total time of potentially decreased photophoretic strength 
is not significant. 

Another source of random rotation to consider is Brownian 
motion. The thermal rotation energy of a spherical particle is 
given by E rot = ^k B T, which leads to 



3k R T„ 



H 1 gas 



45fc B r„ 



./ 



par 



87!T 5 



Ppar 



(31) 



with J par being the sphere's moment of inertia. To rotate by 180° 
and change front and backside it is 



8ttV 5 



Trot(lSO') - 



Ppar 



45k B T : 



(32) 



gas 



For chondrules this is on the order of 10 s. The time scale for 
thermal conduction through the particle is given by 



•cond — i Ppar C par 



(33) 



Here, c par is the particle's heat capacity at constant pressure. For 
chondrules T can d is on the order of 1 s. Therefore, a temperature 
gradient reestablishes itself much faster than thermal rotation 
can disalign it and hence Brownian Rotation can be neglected, 
too. 

Besides random rotation exciting a continuous rotation is 
possible by photophoresis itself. As particles are never ideal 
spheres torques are induced. One can argue though that only 
torques around the direction of light will induce a systematic ro- 
tation around the axis aligned to the direction of light. Any other 
torque wil l only lead to the stable alignm ent of the particle. To 
prove this Ivan Evmeren & Wurml (1201 2l) carried out a detailed 
study of particle rotation. In their setup particles were levitated 
by m e ans of thermophoresis a nd ph otophoresis dKelling et al.l 
1201 lh . Ivan Evmeren & Wurml d2012l) found that at least more 
than 95% of all particles rotate around the vertical direction 
which coincides with the direction of illumination. Such a rota- 
tion does not change front and back with respect to the illumina- 
tion and is not important for the photophoretic strength along the 
direction of light. For this reason rotation is no issue in particle 
transport by stellar radiation. It might be noted that the situation 
is different for radiation pressure acting in the current solar sys- 
tem. However, here no gas and therefore no damping is present. 

Heat exchange with the ambient gas: In our heat transfer 
calculations we did not include heat exchange with the ambient 
gas. Heat transfer to the ambient gas for large Knudsen numbers 
can be quantified as (Sto ffels et alJll996h 



+ 1 



Icond 



aATp. 



2k h 



ft Tgas^gas 



(34) 



8 



C. Loesche and G. Wurm: Thermal and Photophoretic Properties of Dust Mantled Chondrules and Sorting in the Solar Nebula 



where / is the heat flux (power per surface area). 

In the case of molecular hydrogen, which consists of two 
atoms, c p /c v can by approximated by 9/7. AT denotes the tem- 
perature difference between the chondrule and the surround- 
ing gas . For a chond r ule at 2AU in a minimum mass neb- 
ula by lHavashi et al.l (1 19851) we have a radiative flux / = 
341 W/m 2 , p = 0.15 Pa, T gas = 262 K. Maximum temperature 
differences are Ar = 10 K. This leads to I amt i = 35 W/m 2 . This 
is an order of magnitude less than the incoming radiative flux 
and can be neglected in the context of photophoretic sorting of 
chondrules in protoplanetary disks. 

Typical laboratory conditions are / = 20kW/nr (maximum 
Ar = 400 K, see fig.©, p = 1 Pa, m gas = 29 u (air), and 
T gas - 293 K which yields I con d = 315 W/m 2 . This is also 
much lower than the incoming light flux and cooling by gas 
interaction still can be neglected. In case of higher pressures 
(p = 10 Pa . . . 100 Pa) heat exchange with the gas becomes im- 
portant, but at the same time the assumption of photophoresis 
in the free flow regime is also no longer valid. For transitional 
disks assumed to be of low pressure our treatment holds. On 
the other hand higher disk (gas) densities lead to a reduced light 
flux as they are usually assumed to go with higher absorption. 
This also changes the photophoretic strength and transport. In 
the extreme case of an optically thick disk, photophoresis does 
not work. Therefore, high density regions of a disk require an ex- 
tended treatment of photophoresis, but this is beyond the scope 
of this paper. 



7. Conclusion 

In summary there are two options to size sort chondrules with 
dust mant les by photophor esis. The mantle thickness variation 
with size dCarballidoll201 lh naturally induces different effective 
thermal conductivities. Plausible porosity variations with size 
act in the same direction and on the same magnitude. Therefore, 
if dust covered chondrules are subject to a directed radiation they 
likely get size sorted by photophoresis. In the original work by 
W urm & Kraussl (|2006j) the solar nebula was assumed to be op- 
tically thin throughout. This is not likely to be the case for the 
very young nebula. However, an increasing number of transi- 
tion disks is fou nd, which have an optical thin inner clearin g 
up to tens of AU dNaiita et al.ll2007HSiciiia-Aguilar et alJl2008h . 
Some of these show active signs of accretion and therefore the 
op tical thin part is s t ill gas eous. Within such a clearing the idea 
by IWurm & Kraussl d2006l) can naturally be applied. The effec- 
tive thermal conductivities of dust mantled chondrules are 1 to 2 
orders of magnitude lo wer than the thermal con ductivity of bare 
chondrules assumed bv lWurm & Kraussl (120061) . This increases 
the photophoretic strength by a similar factor. Therefore, dust 
mantled chondrules can concentrate in the asteroid belt region 
even if the inner clearing has a significantly reduced gas density. 
Creating the inner clearing takes some time. This would also fit 
well with the fact that chondrules are forming their parent bodies 
only a few milli o n year s after the beginning of the solar nebula. 
lHaack & Wurml (120071) also considered the idea that the inner 
edge of the disk might be important. This would shift the focus 
from the spatially extended inner clearing to the actual edge that 
is moving outwards. Here the disk is supposed to be of its orig- 
inal density. Particles exposed to the starlight at the edge would 
feel a photophoretic force and be pushed outwards until the op- 
tical depth is too large. The final position with respect to the 
edge will depend on the susceptibility to photophoresis and par- 
ticles with lower effective thermal conductivity will move fur- 



ther out. This will introduce a size sorting at the edge of the disk 
as it moves inside out. If parent bodies form in this increased 
density environment chondrules will be size sorted. A detailed 
treatment of this is beyond the scope of this paper as it would 
have to focus on disk models, while this work is dedicated to the 
photophoretic properties needed as input. Here, we evaluated the 
effective thermal conductivity of a dust mantled chondrule in the 
context of photophoretic forces. It strongly depends on the size 
ratio between mantle and core and the thermal conductivity of 
the mantle. Based on these results a size sorting of dust mantle 
chondrules by photophoresis is possible. 
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